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ABSTRACT 

We have tested for luminosity, colour and morphology dependence of the degree of 
filamentarity in seven nearly two dimensional strips from the Sloan Digital Sky Sur- 
vey Data Release Four (SDSS DR4). The analysis is carried out at various levels of 
coarse graining allowing us to address different length-scales. We find that the brighter 
galaxies have a less filamentary distribution than the fainter ones at all levels of coarse 
graining. The distribution of red galaxies and ellipticals shows a higher degree of fil- 
amentarity compared to blue galaxies and spirals respectively at low levels of coarse 
graining. The behaviour is reversed at higher levels of coarse graining. We propose a 
picture where the ellipticals are densely distributed in the vicinity of the nodes where 
the filaments intersect while the spirals are sparsely distributed along the entire extent 
of the filaments. Our findings indicate that the regions with an excess of ellipticals are 
larger than galaxy clusters, protruding into the filaments. We have also compared the 
predictions of a semi-analytic model of galaxy formation (the Millennium Run galaxy 
catalogue) against our results for the SDSS. We find the two to be in agreement for 
the M* galaxies and for the red galaxies, while the model fails to correctly predict the 
filamentarity of the brighter galaxies and the blue galaxies. 

Key words: methods: numerical - galaxies: statistics - cosmology: theory - cosmology: 
large scale structure of universe 



1 INTRODUCTION 

One of the most significant aim of all large redshift sur- 
veys is to determine the spatial distribution of galaxies 
in the Universe. Over the last few decades redshift sur- 
veys have revealed the large scale structures in the Uni- 
verse i n it's full glory (e.g. CfA , iGeller fc Huchralll989t 
LCRS, IShectman etaHll99d 2dFGRS. IColless et al. 1200 ll 
and SDSS, IStoughton et alJ I2002L lAbazaiian et alJ l200l 
lAbazaiian~et^iTl20o3) ~ It is found that the vast major- 
ity of galaxies preferentially reside in an intricate net- 
work of interconnected filaments and walls surrounded by 
voids. The fila ments have a dimension of many tens of 
Mpc in le ngth jBharadwai et al.l |2004 and 1-2 Mpc i n 
thickness (iRatcliffe et alJ I1996L iDoroshkevich et alJl2004) . 
In fact the Local Group in which the Milky Way resides, 
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is believed to lie along an extended filament that origi- 
nates near the Ursa Major and includes both the IC342 
and M81 groups and connects with another primary fil- 
ament which originates fr om the Virgo cluster and in- 
clude s the Sculpture group iPeebles. Phelps. Shava fe Tullvl 
l200lT) . The complex filamentary network of galaxies, the 
so called 'Cosmic web' is possibly the most striking vis- 
ible feature in the large scale structure of the Universe. 
Quantifying the Cosmic web and tracing it's origin has 
remained one of the central issues in cosmology. The 
analysis of filamentary patterns in the galaxy distribu- 
tion has a long histor y dating back to a few papers 
in th e mid-eighties by IZel'dovich , Ei nasto fc Shandarir 
<1982ft . IShandarin fc Zeldovichl dl983Tl and lEinast^^t^rtl 
il984) . More recently, extensive studies of the LCRS re- 
veal a ric h network of interconnec t ed filaments surroundin g 
voids (eg. [ Shandarin fc Yesslli"99gl : iBharadwai et alJl2000) : 
I Miiller et al.ll2000l : IDoroshkevich et alJl200ll : iTrac et all 
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| 2002l : iDoroshkevich et alJ 120041 ; iBharadwai et alJl20o4 - 
lBhara^wa^7^PandB^ ~ l2004ll have compared the filamen- 
tarity in the LCRS galaxy distribution with ACDM 
dark matter N-body simulations to show that the 
two are consistent provided a mild galaxy bias is in- 



cluded. IShetbl <2t)0fl has used a technique SURFGEN 
JSheth et alJ |2003) to study the geometry, topology and 
morpho logy of the superclusters in mock S PSS cata- 
logues. iBasilakos, Plionis, fc Rowan-Ro binsonl (l200lT) and 
lKolokotrCTris'^B^sikIk^s^^PllCTn^l2^02l) have studied the 
super-cluster void network in the PSCz and the Abell/ACO 
cluster catalogue respe ctively, finding filamentarity to be 
the d ominant feature. IPimbblet. Drinkwater fc Hawkriggl 
(2004) have studied the intercluster galaxy filaments in the 
2dFGRS. 

The Sloan Digital Sky Survey (SDSS) jYork et alj2000Tl 
is curr ently the largest galaxy reds hift survey. In a recent 
paper iPandev fc Bharadwaill2005l) we have analyzed the 
filamentarity in the two equatorial strips of the First Data 
Release of the SDSS (SDSS DR1). These strips are nearly 
two dimensional (2-D). We have projected the data onto 
a plane and analyzed the resulting 2-D galaxy distribu- 
tion. The large volume and dense sampling of the SDSS 
allows us to construct volume limited subsamples extending 
over lengthscales which are substantially larger than possible 
with earlier surveys. We find evidence for connectivity and 
filamentarity in excess of that of a random point distribu- 
tion, indicating the existence of an interconnected network 
of filaments. The filamentarity is found to be statistically 
signifi cant up to length-scal es 80/i _1 Mpc and not beyond 
JPandev fc Bharadwall2005ft . Further, we show that the de- 
gree of filamentarity exhibits a luminosity dependence with 
the brighter galaxies having a more concentrated and less 
filamentary distribution. 

It is now quite well accepted that galaxies with 
different physical properties are differently distributed in 
space. Studies over several decades have established that 
ellipticals and spirals are not distributed in the same way. 
The ellipticals are found predominantly inside regular, 
rich clusters whereas the field galaxies are mostly spirals. 
This effect is referred to as "morphological segregation" 
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of the two-point correlation function of galaxies with 
different colours shows the red galaxies, which are mainly 
ellipticals with old stellar populations, to have a stronger 
clustering as comp ared to the blue galaxies which are 
mostly spirals (e.g. | Willm er 1 da Costa, fc Pellegrin 3 Il998t 
iBrown. Webstar, fc Bovld 12001 IZehavi et alJl2005ll . Stud- 
ies of the topology using the genus statistics show that 
red galaxies exhib it a shift towa rds a meatball to pology 
jHovle et alJ l200l SDSS EDR, iPark et alJ l200l SDSS 
DR3) implying that they prefer to inhabit the high density 
environments. It is found that luminous galaxies exhibit 
a stronge r clustering than their fa i nter counterparts (e.g. 
Hamilton! ll99Sl IDavis et alJ Il988t IWhite. Tullv fc Davis I 
19881. IPark. Voeelev. Geller fc Huchral [1994 
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l200ll IZehavi et al-lEool) . The difference inc rease m arkedly 
abov e the char acteristic magnitud e M* llNorberg et alJ 
I2001L 2dFGRS, IZehavi et alTl2005l SPSS DR2) of the 
Schecter luminosity function ( Schcchtcr; ll976fl . The detailed 
luminosity dependence has been difficult to establish 
because of the limited dynamic range of even the largest 

reds hift surveys. 

lEinasto et alJ <l2003f) have studied the luminosity dis- 
tribution of galaxies in high and low density regions of 
the SDSS to show that brighter galaxies are preferentially 
distrib uted in the high density environments. IGoto et alJ 
( 2003) have studied the morphology density relation in the 
SDSS (EDR) and found that this relation is less noticeable 
in the sparsest regions indicating the need for a denser en- 
vironment for fostering the physical me chanisms responsi- 
ble f or the galaxy morpholo gical changes. iHogg et alJ (2003) 
and iBlanton et alJ 11200301 find a strong environment de- 
pendence for both the colour and luminosity for the SDSS 
galaxies. 

The dependence of clustering on galaxy properties 
like the luminosity, colour and morphology provides very 
important inputs for theories of galaxy formation. It is a 
natural prediction of hierarchical structure formation that 
the rarer objects which correspond to peaks in the density 
field s hould be mor e strongly clustered than the populatio n 
itself (iKaiserl 11984. IDavis et alJ Il985l . IWhite et all Il987f) . 
A possible interpretation of the observations is that the 
more luminous galaxies are hosted in higher density peaks 
as compared to the fainter galaxies. In this scenario the 
properties of a galaxy are largely decided by the initial 
conditions at the location where the gal axy is formed. This 
idea i s implemented in the halo model Jing^Mo^^Fiorncr 
I2000L ISeliakl 12000. iPeacock fc Smith 
2001 IScoccimarro fc She 
jBejdnid^^^Veinber^ 
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Scrantonl l200l lYang. Mo fc van den Boschl 12003. . 
Yang. Mo. Jing fc van den Boschl 12003ft where it is postu- 
lated that all galaxies lie in virialised halos and the number 
and type of galaxies in a halo is e ntirely determined by the 
halo mass. lCoorav fc She th (2002) provide a detailed review 
of the halo model. It is well known that galaxy-galaxy 
interactions and environmental effects are also important 
in determining the physical properties of a galaxy. For 
example, ram pressure stripping in galaxy clusters is a 
possible mechanism for the transformation of a spiral 
galaxy into an elliptical. Galaxy-galaxy interactions could 
also provide a mechanism for morphological transformation. 
In the classical picture galaxies evolve in isolation and 
morphology is determined at birth, while in the hierarchical 
model galaxies have no fixed morphology and they can 
evolve depending on the nature of their mergers. A study 
of the spatial distribution of different kind of galaxies holds 
the potential to probe the factors responsible in determining 
galaxy luminosity, colour and morphology. 

Much of the research on the luminosity, morphology 
and colour dependence of the galaxy distribution has fo- 
cused either on the morphology segregation in clusters or 
has studied the behaviour of the two point correlation as a 
function of these galaxy properties. Filaments are the largest 
known statistically significant coherent features. They have 
been shown to be statistically significant to scales as large 
as 80 hT 1 Mpc. In this paper we study the luminosity, colour 
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and morphology dependence of filaments in the Sloan Dig- 
ital Sky Survey Data Release Four (SDSS DR4). This al- 
lows us to study the distribution of different types of galax- 
ies on the largest length-scales possible. Such a study will 
enable us to estimate the length-scale over which the fac- 
tors responsible for the galaxy segregation operate. While 
the two-point correlation function completely characterizes 
the statistical properties of a Gaussian random field, it is 
well known that the galaxy distribution is significantly non- 
Gaussian. In fact the presence of large-scale coherent fea- 
tures like filaments is a clear indication of non-Gaussianity. 
There have been studies of how the thre e point correlation 
function depends on galaxy properties (|Kavo et alJ 120041: 
Ijing fc Bornerll2004l : iGaztanaea et alj l2005). While the first 
two papers do n o t find any statistically significant effects, 
iGaztafiaga et all (120051) find statistically significant colour 
and luminosity dependence. The present work studies how 
the largest known non-Gaussian features, the filaments, de- 
pend on galaxy properties. 

This paper prese nts progress on many counts compared 
to our earlier work iPandev fc Bharadwaill2005T) based on 
the SDSS DR1. The earlier work was restricted to single 
volume limited subsamples from each of the two equatorial 
strips. The absolute magnitude range was divided into two 
equal halves, and these were analyzed to test for a lumi- 
nosity dependence in the filamentarity. For a statistically 
significant conclusion it is necessary to also estimate the 
cosmic variance (sample to sample variation) of the filamen- 
tarity which was done by bootstrap resampling. We found 
that the difference in filamentarity between the two mag- 
nitude bins was in excess of the fluctuations expected from 
cosmic variance, thereby establishing the statistical signifi- 
cance. The SDSS DR4 covers a large contiguous region in the 
Northern Galactic Cap with a few gaps still to be covered. 
The SDSS has been carried out in overlapping strips of 5° 
width oriented along great circles. Since the detection of co- 
herent large-scale features requires a completely sampled re- 
gion with no gaps, we have extracted seven non-overlapping 
strips of width 2° each. The volume corresponding to each 
strip is nearly two dimensional and we have collapsed the 
thickness and analysed the resulting 2-D galaxy distribution. 
The increased number of strips provide a better estimate of 
the cosmic variance leading to more robust conclusions. The 
analysis has been carried out on five overlapping bins in ab- 
solute magnitude. We have extracted separate volume lim- 
ited subsamples corresponding to each magnitude bin and 
compared the filamentarity across them. The earlier analy- 
sis has been extended here to also consider the colour and 
morphology dependence, which has been tested in each of 
the volume limited subsamples. 

Here we would like to mention that the statistics that 
we are using to quantify the filamentarity is not absolute 
in the sense that it is sensitive to the volume and galaxy 
number density of the sample. It is thus only meaningful 
to compare different galaxy samples provided they have the 
same volume (identical shape and size) and number density. 
This is an important consideration when constructing the 
galaxy samples for which we test the luminosity, colour and 
morphology dependence. 

A brief outline of our paper follows. Section 2 describes 
the data and Section 3 the method of analysis. Our results 
for the SDSS data are presented in Section 4, in Section 



5 we compare our results with the predictions of a semi 
analytic model of galaxy formation and finally we present 
our conclusions in Section 6. 

It may be noted that we have used a ACDM cosmologi- 
cal model with fl m o = 0.3, £Iao = 0.7 and h = 1 throughout. 



2 DATA 

The SD SS is an imagin g and spectroscopic survey of 
the sky jYork et al.ll200rjT) in five photometric bandpasses, 
u,g,r,i,z with effective wavelengths of 35 40A°, 476 0A°, 
628QA , 7690A° and 92 50A° respectively dFukugita et alJ 
Il996t ISmith et alJ l2002Tl to a limiting r band magnitude 
~ 22.2 with 95% completenes s. A series of pipe lines that per- 
form astro metric calibration i Pier et al.ll2003r) , photometric 
reductio n iLupto n et al.|l200S ) and photometric calibration 
jHogg et alJl200ll) are used to process the imaging data and 
then objects are selected f rom the imaging data for spec- 
troscopy. In the SDSS DR4 jAdelman-McCarthv et all2005l) 
the imaging data covers 6670 deg 2 of the sky whereas the 
spectroscopy covers a total area of 4783 deg 2 . The spec- 
troscopic data includes 673,280 spectra with approximately 
480,000 galaxies, 64,000 quasars and 89,000 stars. The sur- 
vey area covers a single contiguous region in the Northern 
Galactic Cap and three non-contiguous region in the South- 
ern Galactic Cap. Our analysis is restricted to the Northern 
Galactic Cap region only. 

The samples analyzed here were obtained from the 
SDSS DR4 Skyserver 1 . This is a web interface to the SDSS 
data archive server. The SDSS surveys the sky in overlap- 
ping stripes 2 of width 5° which form great circles on the sky. 
The stripes are most conve niently described usingj survey 
co-ordinates (A, rj) defined in lStoughton et all i2002h . Lines 
of constant r\ are great circles aligned with the stripes and 
lines of constant A are are small circles along the width of 
the stripes. Each stripe is centered along a line of constant r\ 
separated from the adjoining stripe by 2.5°. We have down- 
loaded data from seven overlapping stripes which form a 
nearly contiguous region shown in Figure^ We have selected 
all objects identified as galaxies with extinction corrected r 
band Petrosian magnitude r < 17.77 over the redshift range 
0.01 > z > 0.2. 

The primary data is not contiguous and has a few gaps. 
Since our analysis requires contiguous regions, we have ex- 
tracted seven non-overlapping strips each of width 2° in 
r] and spanning 90° in A (Figure and Table 1). Only 
galaxies with extinction corrected Petrosian r band mag- 
nitude in the range 14.5 > m r > 17.77 were used. For 
each galaxy the absolute magnitude was computed using its 
redshift, a pparent m agnitude , K correction a nd reddening 
correction iSchleee^ Finkbeiner fc DavisHl998l) . We adopt a 
polynomial fit (Park et al. 2005:) for the mean K-correction, 
K(z) = 2.3537(0 - 0.1) 2 + 1.04423(2 - 0.1) - 2.5log(l + 0.1). 

We construct volume limited sub-samples using the ab- 
solute magnitude bins listed in Table 2. The five overlapping 
absolute magnitude bins are referred to as "bin 1" , "bin 2" , . . . 
"bin 5" in order of increasing luminosity. Each bin extends 

1 |http://cas.sdss.org/dr4/en/| 

2 http:/ /www. sdss.org/dr4/coverage/atStripeDef.par 
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Figure 1. The shaded part shows the contiguous region cov- 
ered by the seven overlapping strips constituting our primary 
data. The bounded regions within show the seven 2° wide non- 
overlapping strips used in our analysis. 



over a different redshift range. The thickness of the volume 
corresponding to each bin increases with redshift. We have 
extracted regions of uniform thickness 6/i _1 Mpc. This cor- 
responds to 2° at redshift 0.06, and for some bins it was 
necessary to discard the range z < 0.06 where the thickness 
was less than 6 ft" 1 Mpc. The redshift range which we use 
for final galaxy samples that were analyzed are shown in 
Table 2. The faintest luminosity bin has the smallest area 
(Table 2), and the area of the other bins increases mono- 
tonically with luminosity. A smaller magnitude range was 
chosen for the faintest bin so as to ensure that it covers a 
reasonably large area while maintaining a number density 
comparable to the other bins. The galaxy number density is 
maximum for bin 2 (Table 2) whic h is centered at M* w hich 
has a value M* = -20.44 ± 0.01 l|Blanton et al.ll2003aT l for 
the SDSS. The number density falls for the other bins. The 
number of galaxies for all the strips in each bin are also 
tabulated in Table 1. 

Studies using N-body simulations show that the statis- 
tics that we use to characterize the filamentarity in the 
galaxy distribution is sensitive to the number density as well 
as the shape and size of the volume of the sample. To test 
for luminosity dependence it is necessary to compare differ- 
ent luminosity bins which as seen in Table 2 have different 
number densities and cover different areas. To ensure that 
these factors do not influence our results, we extract regions 
of identical shape and size as bin 1 from the other bins and 
use these to test for luminosity dependence. Tests using dark 
matter ACDM N-body simulations show that galaxy num- 
ber density variations of ~ 50% do not produce statistically 
significant effect on the filamentarity (Figure 0. We have 
discarded randomly chosen galaxies from bins 2 and 3 so as 
to make the galaxy number density of each strip of these 
bins exactly equal to the mean density of bin 1. The galaxy 
number density of bin 4 is within 10% of that of bin 1 while 
it is 27% lower in bin 5. Since these are lower than bin 1, 
and as they are within the permissible range of variation 
we do not dilute the galaxy number density in bins 4 and 5. 
Figure|5|shows the galaxy distribution in two different lumi- 



nosity bins drawn from the same strip. It is to be noted that 
the samples with reduced area and diluted galaxy number 
density are used only for testing luminosity dependence. 

We separately study the colour and morphology depen- 
dence in each of the magnitude bins listed in Table 1 and 
Table 2. We do not compare the colour and morphology de- 
pendence across different luminosity bins. When testing for 
colour dependence in a fixed luminosity bin, all the galaxies 
(Table 1) in the entire area of the bin (Table 2) are classified 
as either red or blue galaxies. The gal axy u — r colours ar e 
known to have a bimodal distribution dStrateva et al .1200 ill . 
In our analysis we determine a value (u — r) c for the colour 
such that it divides each magnitude bin into equal number of 
red (i.e. u~r > (it— r) c ) and blue (i.e. u — r < (u—r) c ) galax- 
ies. As a consequence the number density of red and blue 
galaxies in each luminosity bin are exactly equal and have a 
value half that of the density shown in Table 2. The value of 
(u~r) c varies slightly across the luminosity bi ns, and the val- 
ues ar e listed in Table 2. It may be noted that lStrateva et all 
feOOlf) finds that a color selection criteria (u — r) c = 2.22 
ensures that the red galaxies are 'ellipticals' with 90% com- 
pleteness. Figure [3] shows the distribution of red and blue 
galaxies in one of the subsamples. 

The morphological classification was carried out us- 
ing the concentration index defined as Ci — rgo/rso where 
rgo and 7-50 are the radii containing 90% and 50% of 
the Petrosian flux respectively. This has been found to 
be o ne of the best paramet er to classify galaxy morphol- 
ogy Ishimasak u et aD 1200 J) . Ellipticals are expected to 
have a larger concentration indices than spirals. It was 
found that d ~ 3. 33 for a pure de-Vaucouleurs profile 
jBlanton et alJ l200ll) while d ~ 2.3 for a pure exp onen- 
tial profikT lstrateva et all200ll) . IStrateva et alJ i2001ri show 
that using a cutoff Ci tC = 2.6 divides the sample into ellip- 
ticals and spirals at 83% completeness. For each luminosity 
bin we have chosen a cutoff Ci jC that partitions the galaxies 
into two equal halves, one predominantly ellipticals and the 
other spirals. The values of Ci, c varies across luminosity bins 
and are listed in Table 2. Figure 01 shows the distribution of 
early and late type galaxies in one of the subsamples. 

Finally we note that there is an incompleteness in the 
SDSS survey arising from a restriction which prevents the 
redshift of two galaxies at a very small angular separation 
to be measured. This incompleteness is not expected to in- 
troduce a luminosity, colour or morphology dependence and 
hence we do not take this into account. 



3 METHOD OF ANALYSIS 

Each of the subsamples described in the previous subsection 
have a much greater (16-35 times) linear dimension com- 
pared to their thickness and so can be reasonably treated as 
two dimensional. The 2D galaxy distribution was embedded 
in a 1 /i _1 Mpc x 1 /i _1 Mpc 2D rectangular grid. The galaxy 
distribution was represented as a set of Is and 0s on a 2- 
D rectangular grid by assigning the values 1 and to the 
filled and empty cells respectively. Note that ~ 10% of the 
filled cells have 2 galaxies while the number is substantially 
smaller for 3 galaxies or more, and hence the Is and 0s pro- 
vide a fair representation of the galaxy distribution on the 
length-scales of interest. The grid cells which are beyond the 
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Figure 2. This shows the distribution of faint and bright galaxies in one of the strips in magnitude bin 1 and bin 4 respectively after 
three round of coarse-graining. 
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Figure 3. This shows the distribution of blue and red galaxies in one of the strips in magnitude bin 4 after three round of coarse-graining. 




Mpc/h 

Figure 4. This shows the distribution of early and late type galaxies in one of the strips in magnitude bin 4 after three round of 
coarse-graining. 
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Table 1. This shows the (A, r)) range of the seven non-overlapping strips extracted from the primary data. For each strip the table shows 
the number of galaxies in volume limited subsamples bin 1, bin 2, .. etc. with absolute magnitude and redshift limits given in Table 2. 



Strip number 


A(°) 


v(°) 


bin 1 


bin 2 


bin 3 


bin 4 


bin 5 


1 


-50 < A < 40 


9 < rj < 11 


875 


2016 


2108 


1943 


1623 


2 


-50 < A < 40 


11 < r\ < 13 


839 


1740 


1780 


1695 


1591 


3 


-60 < A < 30 


13 < i) < 15 


663 


1492 


1665 


1597 


1540 


4 


-60 < A < 30 


15 < v < 17 


678 


1547 


1629 


1569 


1593 


5 


-50 < A < 40 


21.5 < r] < 23.5 


846 


1878 


1863 


1703 


1595 


6 


-50 < A < 40 


24 < T) < 26 


782 


1905 


1854 


1839 


1668 


7 


-50 < A < 40 


26 < r) < 28 


648 


1550 


1573 


1692 


1633 



Table 2. This shows the absolute magnitude and redshift limits for the different volume limited subsamples analyzed. The first three 
bins actually extend to lower redshifts, but only the region beyond z = 0.06 has been kept so that it is possible to extract a volume of 
uniform thickness 6h~ 1 Mpc. The area and the average galaxy number density with 1 — cr variations from the 7 strips are also shown. 
The last two column shows the value of the galaxy colour (u — r) c and the value of the concentration index Ci jC used to divide the data 
into equal numbers of red/blue and elliptical/spiral galaxies respectively. 



bin 


Absolute Magnitude range 


Redshift range 


Area [ 10 4 /i" 2 Mpc 2 ] 


Density [10 -2 /i 2 Mpc~ 2 ] 


(u - r) c c i>c 


bin 1 


-19.5 > M r > -20 


0.06 < z < 0.093 


3.37 


2.25 ± 0.26 


2.14 2.53 


bin 2 


-19.75 > M r > -21.25 


0.06 < z < 0.103 


4.66 


3.71 ± 0.41 


2.35 2.67 


bin 3 


-20 > M r > -21.5 


0.06 < z < 0.114 


6.22 


2.86 ± 0.27 


2.41 2.72 


bin 4 


-20.25 > M r > -21.75 


0.06 < z < 0.126 


8.08 


2.13 ± 0.15 


2.45 2.75 


bin 5 


-20.5 > M r > -22 


0.067 < z < 0.14 


9.87 


1.63 ±0.037 


2.5 2.78 



boundaries of the survey were assigned a negative value in 
order to distinguish them from the empty cells within the 
survey area. 

The next step is to use an objective criteria to identify 
the coherent large-scale structures visible in the galaxy dis- 
tribution. We use a "friends-of-friends" (FOF) algorithm to 
identify interconnected regions of filled cells which we refer 
to as clusters. In this algorithm any two adjacent filled cells 
are referred to as friends. Clusters are defined through the 
stipulation that any friend of my friend is my friend. The dis- 
tribution of Is on the grid is very sparse with only ~ 1% of 
the cells being filled. Also, the filled cells are mostly isolated, 
and the clusters identified using FOF, which contain only a 
few cells each, do not resemble the large-scale coherent struc- 
tures seen in the SDSS strips. It is necessary to coarse-grain 
the galaxy distribution so that the large scale structures may 
be objectively identified. In every iteration of coarse-graining 
we fill up all the empty cells adjacent to every filled cells (i.e. 
cells at the 4 sides and 4 corners of a filled cell) , causing every 
filled cell to grow fatter. The size of an isolated filled cell is 
(2A^ + 1) x (2N + 1) after N iterations of coarse-graining and 
for example it is 7/i _1 Mpc x 7/i _1 Mpc after 3 iterations. 
This causes clusters to grow, first because of the growth of 
filled cells, and then by the merger of adjacent clusters as 
they overlap. Increasing the coarse-graining further induces 
the percolation transition when a large fraction of the clus- 
ters connect up into a single, multiply connected network of 
filaments encircling voids. As coarse- graining proceeds even 
further the network grows to some extent, and eventually 
the elements become very thick filling up the entire survey 
region washing away any visible pattern. The filling factor 
FF, defined as the fraction of cells within the survey area 
that are filled, le. 



FF 



Total No. of Filled Cells 



increases from FF ~ 0.01 to FF = 1 as the coarse-graining 
is increased. The filling factor is around FF ~ 0.3 after 3 
rounds of coarse- graining (Figure|HJl. So as to not restrict our 
analysis to an arbitrarily chosen level of coarse-graining, we 
analyze the clusters identified in the pattern of Is and 0s 
after each iterat ion of coarse-graining. Earlie r studies using 
the SDSS DR1 JPandev fc Bharadwaill2005i) show that the 
percolation transition occurs in the FF range 0.5 — 0.6. 

The geometry and topology of a two dimensional cluster 
can be described by the three Minkowski functionals, namely 
its area S, perimeter P, and number of holes or genus G. 
We have tested that the cluster area S is proportional to 
the actual number of galaxies within the cluster boundary. 
This holds at each level of coarse-graining, with the propor- 
tionality constant increasing with coarse-graining. The ratio 
T = S/P characterizes the thickness of the cluster, and P 
the extent of its boundary. Since the clusters are largely fil- 
amentary (as we shall see later), we may interpret P as an 
estimate of the length of the cluster. The ratio P/G has 
dimensions of length an d it is particularly sign ificant after 
the onset of percolation (Bh aradwai et alH 2000'l as it char- 
acterizes the length of the cluster boundary per void. It is 
possible to quantify the shape of the cluster using a single 
2D "Shapefinder" statistic (Bharadwaj et al. 2000) which is 
defined as the dimensionless ratio 



T = 



47rS' 



4715" 



(2) 



Total No. of Cells Inside the Survey Area 



(1) 



which by construction has values in the range < T < 
1. It can be verified that T = 1 for an ideal filament which 
has a finite length and zero width, whereby it subtends no 
area (S = 0) but has a finite perimeter (P > 0). It can be 
further checked that T = for a circular disk, and interme- 
diate values of J- quantifies the degree of filamentarity with 



Luminosity, colour and morphology dependence of galaxy filaments in SDSS DR4 7 



the value increasing as a cluster is deformed from a circular 
disk to a thin filament. 

The definition of T needs to be modified when working 
on a rectangular grid of spacing I. An ideal filament, repre- 
sented on a grid, has the minimum possible width i.e. I, and 
its perimeter P and area S are related as 2S = (P — 21)1. At 
the other extreme we have P 2 = 165" for a square shaped 
cluster on the grid. We introduce the 2D Shapefinder statis- 
tic 



T ■ 



(P 2 - 165) 
(P - 4Z) 2 



(3) 



to quantify the shape of clusters on a grid. By definition 
0< T < \. T quantifies the degree of filamentarity of the 
cluster, with T = 1 indicating a filament and T = 0, a 
square, and T changes from to 1 as a square is deformed 
to a filament. 

At every stage of coarse-graining we have a number of 
clusters, each of which is characterized by its shape T% and 
size Si. It is desirable to identify a single statistical quantity 
to capture the filamentarity of the entire set of clusters. We 
have considered two possibilities Fi and F2, respectively the 
first and second area weighted moments of the filamentarity. 
Our aim being to quantitatively compare the filamentarity 
of different galaxy samples and draw statistically significant 
conclusions, we find that it is more advantageous to use F2 
as a statistical discriminator. Though both Fi and F% show 
similar behaviour, the differences in F2 between point dis- 
tributions with different clust ering properties is mor e pro- 
nounced as compared to Fi jBharadwai et al]l2000li . The 
average filamentarity F2 is defined as the mean filamentar- 
ity of all the clusters in a slice weighted by the square of the 
area of each clusters 



F 2 



(4) 



In the current analysis, we study the average filamentar- 
ity F2 as a function of FF to quantify the degree of filamen- 
tarity in each of the SDSS subsamples and investigate how 
the average filamentarity F2 changes with intrinsic galaxy 
properties like luminosity, color and morphology. We also 
use the area weighted second moment (defined analogous to 
eq. (QJ) of S/P, P and P/G to quantify the average cluster 
thickness, length and length per void respectively. 

At a specified level of coarse-graining the value of the 
filling factor (FF) has small differences from sample to sam- 
ple as seen in Figure^] This makes it difficult to compare the 
filamentarity in different samples using the values of the av- 
erage filamentarity (F2) as a function of FF (Figure |SJ. We 
overcome this by interpolating the values of F2 as a function 
of FF at uniformly chosen values which are same for all the 
samples. The average thickness and length are also interpo- 
lated and analyzed in a similar fashion. Further, since the 
later stages of coarse graining tend to wash out the features 
of the galaxy distribution, we have restricted our analysis to 
the range of filling factor < FF < 0.75. 



4 RESULTS 

We first study how the average filamentarity varies across 
the different luminosity bins discussed in Section 2. For each 




3 4 5 6 

Number of coorse-qroininq 

Figure 5. This shows the Filling Factor (FF) as a function of 
the number of iterations of coarse-graining applied to the initial 
galaxy distribution for seven SDSS different strips in the lumi- 
nosity bin indicated in the figure. 




0.3 0.4 0.5 
Fillinq Factor 

Figure 6. This shows the average filamentarity (F2) as a func- 
tion of filling factors (FF) for all seven individual SDSS strips in 
the luminosity bin indicated in the figure. The markers show the 
values for each iteration of coarse-graining, whereas the curves 
show the interpolation used to determine F2 at equally spaced 
values of FF. 



luminosity bin we have seven realisations of the galaxy dis- 
tributions from the seven non-overlapping strips. The results 
from all the individual strips in luminosity bin 1 are shown 
together in Figure |S| For each bin we use the results from 
these seven samples to determine the mean and variance of 
F2 at uniformly chosen values of FF. The results are shown 
in the left panel of Figure |7| We find that for a fixed value 
of FF, the value of F2 decreases with increasing luminosity. 
We note that there is some deviation from this behaviour at 
FF > 0.4 between bins 3 and 4, but the overall trend is still 
valid when we compare these bins with fainter or brighter 
samples. Further, the luminosity dependence is enhanced 
with increasing luminosity. 

At each values of FF we use the Student's t-test to 
quantify the statistical significance of the difference in the 
mean F2 between different luminosity bins. The variance of 
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Figure 7. The average filamentarity (.Fa) as a function of filling factor (FF) for the galaxies divided by luminosity (left), colour (middle) 
and morphology (right). The 1 — a errorbars are shown only for a single curve in each panel, the other curve(s) shown in the same panel 
have similar errorbars. The different magnitude bins used in studying the luminosity dependence (Table 2.) are shown in the left panel. 
The colour and morphology dependence are shown for a single luminosity bin (bin 4 of Table 2.), the other luminosity bins show similar 
colour and morphology dependence. 



F2 is very similar in the different luminosity bins and we 
estimate the standard error for the difference in the means 
using 

fy^ (xi~x A ) 2 +y^ (xi-x B ) 2 " " 

D y N A +N B -2 \N A ^N B > 

where the sum is over the points in the two samples A 
and B which are being compared. Here xa and xb referes 
to the mean, and Na and Nb refers to the number of data 
points. In our case JVa = N B = 7 ■ We use t = ^±b_ 

so 

to estimate the significance of the differences in the means. 
This is expected to follow a Student's t-distribution with 12 
degrees of freedom. We accept the difference in the means 
as being statistically significant if the probability of its oc- 
curring by chance is less than 5%. We find that there are 
no statistically significant differences between adjacent lu- 
minosity bins except the two brightest ones (bins 4 and 5). 
All the non-adjacent bins (eg. bin 1 and bin 3, etc.) show 
statistically significant differences at most values of FF. As 
noted earlier, the F2 curves are sensitive to the galaxy num- 
ber density. While bins 2 and 3 have been culled so that they 
have the same number density as bin 1, bins 4 and 5 have 
number densities which are 10% and 27% lower respectively. 
Figure [S] shows the effect of varying the number density by 
50% in 18 independent realizations of mock galaxy samples 
identical to bin 1 in thickness, area and number density (Ta- 
ble 2) drawn from dark mat ter ACDM N-body simulations 
(Bharadw aj fc Pandevl20o4) . The test shows that 50% vari- 
ations in the number density do not introduce a statistically 
significant effect on the filamentarity except at FF = 0.05. 
This is quite distinct from the luminosity dependence seen 
in the left panel of Figure [7| where there are statistically sig- 
nificant differences among the different luminosity bins at 
all FF except for FF = 0.05. This clearly establishes the 
luminosity dependence to be genuine and not an artifact due 
to number density variations. 
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Figure 8. This shows how the average filamentarity (Fa) as a 
function of filling factors (FF) changes if the galaxy number den- 
sity is changed by 50%. Each curve shows results from 18 strips 
identical in area and thickness as bin 1 (Table 2) drawn from 
3 independent ACDM dark matter N-body simulation. The solid 
curves and the errorbars are for mock strips whose galaxy number 
density exactly matches the mean density of bin 1 (Table 2). 

We have separately analyzed the filamentarity of the 
red and blue galaxies in each luminosity bin. The results 
for bin 4 are shown in the middle panel of Figure Q We find 
that there is a statistically significant colour dependence. At 
the smallest value of FF the red galaxies exhibit a higher 
filamentarity than the blue ones. There is a cross-over at 
FF ~ 0.2, and the blue galaxies have a higher filamentarity 
than the red ones at larger values of FF. A similar trend is 
found in all the other luminosity bins, results for which are 
not shown here. 

The morphology dependence has been studied sepa- 
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rately for each luminosity bin. The results for bin 4 are 
shown in the right panel of Figure |7| The ellipticals have 
a higher F2 as compared to spirals for FF < 0.25 whereas 
the spirals have a higher F2 at larger values of FF. We find 
that the differences are statistically significant for FF < 0.1 
and FF > 0.4. The results are similar for the other lumi- 
nosity bins not shown here. 

Figure [5] shows the average cluster length, length per 
void and thickness as a function of FF for two different lumi- 
nosity bins, and also the colour and morphology dependence 
for a fixed luminosity bin (bin 4). We see that the clusters are 
longer for the faint galaxies, blues galaxies and the spirals as 
compared to the bright galaxies, red galaxies and ellipticals 
respectively. Unlike the length which increases nearly mono- 
tonically with coarse- graining, the cluster length per void is 
found to be quite stable to coarse- graining after the onset of 
percolation. The luminosity, colour and morphology depen- 
dence found in the cluster length is exactly reversed when 
the length per void is considered. We also note that the lu- 
minosity, colour and morphology dependence of the length 
and the length per void are not at a high level of statisti- 
cal significance, except at a few points where the data are 
seen to lie outside the 1 — a errorbars. The average cluster 
thickness, on the other hand, shows statistically significant 
differences at all values of the filling factor. At the same 
value of FF, the clusters are thicker for the bright galax- 
ies, red galaxies and elliptical galaxies as compared to their 
faint, blue and spiral counterparts respectively. 



5 COMPARISON WITH A SEMI ANALYTIC 
MODEL OF GALAXY FORMATION 

Semi analytic models have, over the last two 
decades emerged as a very powerful to ol for 
studying galaxy formation (| Wlii^^^F^nk] ^9j|- 
Lacev fc Silklll99ll iKauffmann. White fc Guiderdonilll99 



Kauffniann fc White! Il993l iLacev et all Il993t ICole et alJ 
1994, Kauffniann 1996, Kauffnia nn. Nusser fc S teinmctz 
~ Tl99a ISomerville fc Primackl [ 1999^ 
ll99l^7Kauffmann et al.1 Il999 



iKai 
jBa 



BauEh^^akl 

Somerville fc Kolat tj 



Cole et alJl200Ci iBenson et al-lEool ISpringel et alJ 12005ft . 

These models describe the evolution of galaxies in a 
hierarchical clustering scenario incorporating all relevant 
physics of galaxy formation processes (eg. gas cooling, 
star formation, supernovea feedback, metal enrichment, 
merging etc) often in an approximate and adhoc fashion. 
The detailed physics of star formation and it's regulation 
by different feedback mechanisms is still poorly understood. 
These models serve as simplified simulations of the galaxy 
formation processes. The output of the semi analytic mod- 
els are statistical predictions of galaxy properties at some 
epoch and the precision of these predictions are directly 
related to the accuracy of the input physics. In this Section 
we investigate the luminosity, and colour dependence of 
galaxy filaments in a particular semi-analytic model and 
compare the findings with those from the SDSS. We have 
used the semi an alytic galaxy catalogs from the Millennium 
Run simulation jSpringel et al.l 12005^1 , one of the largest 
simulation of the growth and evolution of cosmic structures 
in the universe. The det ails of the simulation can be found 
in ISpringel et alJ (|2005). The physical treatment of the 



galaxy forma t ion pr ocesses in this model are described in 
ICroton et alJ ll200fj) . The spectra and magnitude of the 
model ga laxies were computed usi ng population synthesis 
models of iBruzual fc Cha riot ( 2003.) and we use the catalog 
where the galaxy magnitudes are available in SDSS u,g,r,i,z 
filters. The catalog contains about 9 million galaxies in the 
full simulation box. We use these to construct mock strips 
from the simulation which are otherwise identical to the 
SDSS strips which were analyzed. 

The luminosity dependence in the semi-analytic model 
is shown in the left panel of Figure 1101 This is found to 
exhibit a strong luminosity dependence. Although the be- 
haviour is qualitatively similar to the actual data, the lu- 
minosity dependence is noticeably less pronounced for the 
brighter galaxies in the semi-analytic model. The differences 
between the behaviour in the actual data and the semi- 
analytic model are clearly seen in the middle panel of Figure 
1101 which, for two luminosity bins, shows the results for the 
actual data and the semi-analytic model together. Note that 
while there is a reasonable agreement for the faint luminos- 
ity bins, there are sharp differences between the data and 
the simulations in the filamentarity of the high luminosity 
galaxies. 

The right panel of Figure 1101 shows the colour depen- 
dence of the filamentarity in the semi-analytic model. Here 
again, although the model is qualitatively similar to the 
data, the colour dependence predicted in the model is sub- 
stantially in excess of those seen in the actual data. The 
differences between the data and the simulation is partic- 
ularly noticeable at low FF where the filamentarity of the 
blue galaxies in the model is found to be substantially below 
that of the actual data. The simulation results for the red 
galaxies are in rough agreement with the SDSS data. 

The concentration index (a) was not directly available 
in the particular simulated catalogue which we have used 
and hence we did not study the morphology dependence for 
the semi-analytic model. 



6 SUMMARY AND CONCLUSION 

Filaments are the largest known statistically significant co- 
herent features visible in the galaxy distribution. We test if 
the degree of filamentarity depends on the galaxy proper- 
ties. We find evidence for statistically significant luminosity, 
colour and morphology dependence in the average filamen- 
tarity F2 studied as a function of the filling factor FF. 

Comparing the average filamentarity F2 of different lu- 
minosity bins (Figure [7J , we find that the fainter galaxies 
have a more filamentary distribution as compared to the 
brighter ones. These differences exist at nearly all values of 
FF. The drop in filamentarity with increasing luminosity is 
found to be particularly more pronounced in the two bright- 
est luminosity bins which we have considered, both of which 
are brighter than M* . Comparing the thickness and length 
(Figure 03, we find that at a fixed value of FF structures 
traced by the bright galaxies are thicker and shorter than 
those traced by the faint ones. While the differences in the 
length are statistically significant only at large FF, the dif- 
ferences in the thickness are statistically significant at all 
values of the filling factor. With increasing coarse-graining 
there is a percolation transition at FF ~ 0.5 beyond which 
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Figure 9. The upper panels show the average cluster length and length per void while the lower panels show the average cluster thickness. 
The length increases nearly monotonically with FF whereas the length per void is stable after the onset of percolation at FF ~ 0.5. The 
length per void values have been divided by factor of 2 for convenience of plotting as the curves overlap otherwise. The left panels show 
the brightest and faintest luminosity bins, while the central and right panels show the colour and morphology dependence respectively 
for a fixed luminosity bin (bin 4). Note that in each panel we have shown 1 — cr error bars for one of the curves only, the other curve has 
similar errorbars. 



there is an interconnected network of filaments encircling 
voids. After the onset of percolation the length per void gives 
an estimate of the circumference of the filaments encircling 
the voids and dividing this by ir gives an estimate of the av- 
erage void diameter. We find that the faint galaxy distribu- 
tion is more porous (larger number of holes or voids) and has 
smaller voids with a typical void diameter ~ 80±62 Mpc 
as compared to the bright galaxy distribution which has a 
typical void diameter ~ 165 ± 140 ft -1 Mpc. Note that these 
estimates of the void diameter could be somewhat flawed 
because the estimated void diameters exceed the radial ex- 
tent (96 h~ x Mpc) of the samples which we have used to 
test luminosity dependence (Figure but they serve the 
purpose of demonstrating the luminosity dependence of the 
void size. Our findings indicate that galaxies of different 
luminosity are not uniformly distributed along the cosmic 
web. The brighter galaxies are preferentially distributed in 
more compact and thicker regions with large voids in be- 
tween whereas the fainter galaxies inhabit thin, elongated 
regions with more numerous voids of smaller diameter. We 
note that stu dies using the projected two-point corre lation 
function (eg. iNorberg et"aHl2002t IZehavi et afll2005h show 
an increase in correlation amplitude with increasing lumi- 



nosity, the effect being pronounced beyond L* . These earlier 
results are consistent with our findings but they do not have 
any information ab out the shapes of the clustering patt erns. 
Recent studies (eg. ICoil et alJkoOSt iPollo et al.| 120051) also 
show evidence for luminosity dependence at higher redshifts 
(*~1). 

The results for the colour and morphology dependence 
of the average filamentarity are somewhat different from 
those for the luminosity dependence. We find that the red 
galaxies and the ellipticals have a higher filamentarity as 
compared to the blue galaxies and spirals respectively at low 
filling factors (FigurcQIJ. There is a cross-over at FF ~ 0.15 
for the colour dependence and FF ~ 0.25 for the morphol- 
ogy dependence. The blue galaxies and the spirals have a 
higher filamentarity as compared to the red galaxies and 
the ellipticals respectively at large filling factors. Consider- 
ing the length and the thickness (Figure [5J , we find that 
the length shows a statistically significant colour and mor- 
phology dependence only at large FF (beyond percolation) 
where the distribution of blue galaxies and the spirals has 
a greater length as compared to red galaxies and ellipticals 
respectively. The red galaxies and ellipticals have a thicker 
distribution, the differences being statistically significant at 



Luminosity, colour and morphology dependence of galaxy filaments in SDSS DR4 11 




0.2 0.4 0.6 0.2 0.4 0.6 0.2 0.4 0.6 0.8 

Filling Factor (FF) 



Figure 10. The average filamcntarity (-F2) as a function of the filling factor (FF) for the galaxies in the Semi-Analytic Model (Millennium 
Run). The left panel shows the luminosity dependence. The middle panel shows the faintest and brightest luminosity bins for both the 
actual data (SDSS) and the Semi- Analytic Model (SAM). The right panel shows the colour dependence for a single luminosity bin (bin 
4 of Table 2.). The 1 — a errorbars are shown only for a single curve in each panel, the other curve(s) in the same panel have similar 
errorbars. 



all FF. Considering the voids, we find that the distribu- 
tion of red galaxies has fewer voids which are larger (di- 
ameter ~ 127 ± 46 h" 1 Mpc) whereas the distribution of 
blue galaxies has more voids which are smaller (diameter 
~ 98 ± 17 /i -1 Mpc). The ellipticals and spirals show simi- 
lar differences with void diameters ~ 115 ± 50 ft -1 Mpc and 
~ 90 ± 21 hT x Mpc respectively. These estimates of the void 
diameter are quite reliable as the sample size (bin 4) is larger 
than the void diameter (Figure [3] and 2J. 

We note that our estimates are somewhat larger 
compared to results using other methods. For example, 
lEl-Ad fc Piraiil <ll997l) have found that the voids have a scale 
of around 40 h' 1 Mpc in IR AS 1.2 Jy redsh ift sur vey. A simi- 
lar conclusion is reached bv lLindner et alJ Q1995) from their 
analysis of the Northern Local Void. A r e cent a nalysis of 
voids in the 2dFGRS by iHovle fc Vogelevl J2004D find that 
voids have typical diameters of ~ 30 hT 1 Mpc. Most of these 
estimates use the area or volume to determine the void di- 
ameter whereas we use the void perimeter. Substructures 
in the void perimeter is possibly the dominant reason why 
we get a larger estimate of the diameter. Further, the large 
variance in the void diameter seems to indicate that they 
have a large spread in sizes. 

The colour and morphology are very strongly correlated 
galaxy properties, the red galaxies being predominantly el- 
lipticals and the blue ones spirals. It is well known that 
ellipt ical are fo und primarily in dense groups and clusters 
l|Dresslerl fl980) whereas the spirals are distributed in the 
field. Our finding that at large length-scales the structures 
traced by the ellipticals are thicker, and less filamentary as 
compared to the spirals is consistent with a picture where 
the entire galaxy population is distributed along an inter- 
connected network of filaments encircling voids, the Cosmic 
Web. The ellipticals preferentially inhabit the dense groups 
and clusters which can be identified with the nodes of the 
Cosmic Web, the places where filaments intersect. The spi- 



rals, on the other hand, are sparsely distributed along the 
filaments. The colour and morphology are also known to be 
correlated with the luminosity. The more luminous galaxies 
are predominantly ellitpicals and the fainter ones spirals. 

The fact that at small filling factors (which we may as- 
sociate with small length-scales) the eliipticals have a more 
filamentary distribution as compared to spirals, whereas the 
opposite is found at large length-scales is quite intriguing. 
We discuss below a possible interpretation of this finding. 
The higher filamcntarity of the ellipticals at low FF seem 
to indicate that in addition to residing in the clusters at the 
intersection of filaments, the ellipticals also extend a little 
along the filaments which originate from the clusters. The 
spirals on the other hand, are distributed along the entire 
extent of the filaments. This is shown schematically in Fig- 
ure ^2 In this picture the ellipticals have a more compact 
distribution, and as a consequence they connect up to form 
filamentary clusters at the initial stages of coarse graining. 
Since these structures are localized near the nodes of the 
filaments and they do not extend along the entire length 
of the filaments, their filamentarity increases slowly during 
the later stages of coarse graining. The spirals, on the other 
hand, are sparsely distributed along the entire length of the 
filaments. They connect up only at a later stage of coarse 
graining to define the entire filamentary network of the Cos- 
mic Web. This is a possible explanation why at low filling 
factors the ellipticals have a higher filamentarity than the 
spirals. Further, it also indicates that with increasing filling 
factor (later stages of coarse graining) the filamentarity of 
the spirals should grow faster than that of the ellipticals as 
seen in (Figure 0. Finally we note that the qualitative pic- 
ture presented here is a plausible model which is consistent 
with the quantitative findings of this paper. We do not claim 
that it is the only possible explanation, and it is presented 
here more in the spirit of a hypothesis rather than a con- 
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Figure 11. This shows our picture for the relative distribution 
of ellipticals and spirals along the Cosmic Web. The ellipticals 
densely inhabit the nodes and are distributed within the encircled 
regions shown in the figure. The spirals are sparsely distributed 
along the larger volume of the entire filaments. 



elusion. Further work is required to establish or refute this 
picture. 

As already mentioned several times, it should be noted, 
that the curves showing F 2 as a function of FF are not ab- 
solute. They depend on the geometry of the volume (both 
shape and size) and the galaxy number density. It is only 
meaningful to compare different galaxy samples for which 
these quantities are the same. Once these are fixed, we can 
attribute changes in the filamentarity to factors like lumi- 
nosity, colour or morphology which we are testing for. 

Observational findings indicate a c lose connectio n be- 
tween the local density and galaxy type jDresslerll980l) with 
an increase in the elliptical and SO population and a decrease 
in the spirals in galaxy clusters. Our findings are consistent 
with this, and it indicates that there is an excess of ellip- 
ticals on scales larger than the clusters extending to some 
extent into the filaments. In the Zel'dovich pancake scenario 
filaments form at the intersection of sheets and clusters form 
at the nodes of the filaments. Structure formation occurs in 
a hierarchical fashion with sheets forming first, matter flows 
along the sheets into the filaments which in turn drain into 
the clusters. The higher galaxy density and a denser hot 
gas environment are possibly the factors leading to a higher 
elliptical fraction in the vicinity of the nodes. 

The relation between our findings and the various 
theories for galaxy formation is an important issue. We 
have compared our findings against the Millennium Run 
JSprineel et al.ll2005l) which is one of the largest simulated 
galaxy catalogues that incorporates the theory of galaxy for- 
mation in a semi-analytic fashion. The semi-analytic model 
has a number of parameters which have been tuned to match 
various observed properties of the local galaxy distribution 
like the luminosity functions and the Tully-Fisher relation. 
In this paper we have compared the filamentarity of the 



galaxy distribution in the semi-analytic model against the 
actual data from the SDSS DR4. Studying this in different 
luminosity bins we find that the two are consistent for L* 
galaxies, while the brighter galaxies show significant differ- 
ences. For the brighter galaxies, at large filling factors the 
semi-analytic model predicts a substantially higher filamen- 
tarity as compared to the actual data, whereas it under- 
predicts the same at small FF. Though the semi-analytic 
model correctly predicts the number density of very lumi- 
nous galaxies as reflected in the luminosity function, it fails 
to correctly reproduce the geometry of their spatial distri- 
bution. For these galaxies the model does not predict a 
sufficiently filamentary distribution at small length-scales 
whereas there is an excess filamentarity predicted at large 
scales. 

Comparing the filamentarity predicted by the model 
against the actual data for galaxies of different colours , 
we find that the two are consistent for red galaxies whereas 
the distribution of the blue galaxies are quite different. For 
these galaxies the model markedly underpredicts the fila- 
mentarity at small length-scales or filling factors. We spec- 
ulate that the discrepancy may be a consequence of two 
possible ingredients of the model, the first being a radio 
mode feedback from AGNs in massive halos introduced to 
stop cooling flows. This effectively quenches star formation 
in the galaxies in these halos thereby transforming them 
into red galaxies. The second p ossibility is the prescript ion 
(plane parallel slab model) of iKauffmann et al.l lll999t) to 
include the effects of dust when calculating the galaxy lu- 
minos ities and colours. It may be noted that ISpringel et ail 
(2005) have found discrepancies in the colour dependence of 
the two point correlation function. They find that the dif- 
ferences in the correlation amplitudes of the red and blue 
galaxies predicted by the model are in excess of that seen in 
the 2dFGRS and SDSS. 

In conclusion we note that the filamentary pattern seen 
in the galaxy distribution exhibits statistically significant 
luminosity, colour and morphology dependence. We have 
speculated on a geometrical picture which can explain some 
features of the colour and morphology dependence of the fil- 
amentarity. Explaining the origin of the observed luminosity, 
colour and morphology dependence in terms of the theories 
of galaxy formation is an important issue which needs to be 
addressed in the future. 
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